Manifestation of multiband optical properties of MgB2 
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The optical conductivity of MgB2 has been determined on a dense poly crystalline sample in the 
spectral range 6 meV - 4.6 eV using a combination of ellipsometric and normal incidence reflectivity 
measurements. ai{u!) features a narrow Drude peak with anomalously small plasma frequency (1.4 
eV) and a very broad "dome" structure, which comprises the bulk of the low-energy spectral weight. 
This fact can be reconciled with the results of band structure calculations by assuming that charge 
carriers from the 2D a-bands and the 3D 7r-bands have principally different impurity scattering rates 
and negligible interband scattering. This also explains a surprisingly small correlation between the 
defect concentration and Tc, expected for a two-gap superconductor. The large 3D carrier scattering 
rate suggests their proximity to the localization limit. 



PACS numbers: 74. 70. Ad, 78.20.Ci, 78.30.-j 

The recently discovered copperless superconductor 
MgB2 ^ has attracted substantial attention. Most re- 
searchers agree that MgB2 is a conventional sp-metal, 
where a combination of strong bonding and sizable den- 
sity of states produces a high critical temperature within 
the standard theory [^|-Q . There is an agreement among 
theorists that there are two qualitatively different sys- 
tems of bands in MgB2: quasi-2D cr-bands, and 3D tt- 
bands, and that the former strongly couple to optical 
E2g phonon at « 600 cm~^. There is also experimental 
pi and theoretical S evidence of two different supercon- 
ducting gaps in this compound. 

Optical spectroscopy is a direct probe of the electronic 
structure, and it has not been supportive of the con- 
ventional picture so far: the reported plasma frequency 
0g] is w 1.5 eV, in utter disagreement with LDA cal- 
culations of ~ 7 eV [§,|[||. It is tempting to ascribe 
this to correlation effects, so that the missing spectral 
weight is shifted to some incoherent excitations. How- 
ever, there is no reason to expect strong correlations in 
MgB2, which does not have any d- or f- elements, and 
is in no proximity to magnetism. Furthermore, recent 
angle-resolved photoemission data on a single crystal Q 
have found a remarkable correspondence of the observed 
bands to LDA calculations, implying that there is virtu- 
ally no renormalization of Fermi velocity in this system, 
as opposed, for instance, to high-Tc cuprates. Different 
values for the carrier scattering rate at T ^ Tc have been 
reported from transmission spectroscopy on thin films: 
150 cm-i (18.6 meV) M, 300 cm"! (37.2 meV) § and 
700 cm~^ (86.8 meV) Q. Lacking large single crystals, 
a study of polycrystalline samples is justified as a first 
step to understand optical properties of this material. In 
this Communication we report the optical conductivity 



of MgB2 obtained from a dense polycrystalline sample 
from ellipsometry and normal-incidence reflectivity in a 
range from 6 meV to 4.5 eV, which comprises both in- 
traband and low-energy interband excitations. We also 
report theoretical calculation of the same quantity, in- 
cluding the electron-phonon interaction (EPI) and the 
polycrystalline nature of the sample. We suggest that 
the intraband conductivity consists of two qualitatively 
different components: a narrow Drude peak (DP) due to 
the (T-bands, and an overdamped DP due to the 7r-band 
electrons. The latter scatter so strongly in our sample 
that one may expect deviations from the Drude shape due 
to localization effects, and indeed we find that the fea- 
ture in question deviates from the Drude formula. This 
model explains naturally a yet unresolved paradox of the 
two-band scenario [g|,|2| : the absence of a pair- breaking 
by paramagnetic impurities, predicted for two-gap super- 
conductors |13[| . In this Communication we concentrate 
on the broad-range optical properties keeping the FIR 
study of the superconducting gap for a separate publica- 
tion. 
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FIG. 1. The DC resistivity (a), the magnetic susceptibility 
(b) and the X-ray difTraction (c) of the MgB2 sample. 

The sample was synthesized from high-purity pow- 
dered Mg and B, as described elsewhere ffl, and then 
treated in a high-pressure BN cell at 1623 K and 5.5 
GPa for 30 minutes and finally quenched to room temper- 
ature. A very dense metallic-like sample was obtained. 
The Tc is about 39 K with transition width of less than 
1 K both from the DC resistivity and the magnetic sus- 
ceptibility (FigH (a, b)). The estimated superconduct- 
ing volume fraction is about 60 %. The X-ray diffrac- 
tometry (Fig.|l|(c)) revealed presence of impurity phases 
which is difficult to avoid in a bulk polycrystalline sam- 
ple. One can estimate that the sample contains some 
5 - 7 % of MgO and about 1 % of other phases. The 
R(300K)/R(50K) ratio is 2.2. 

The sample was cut and dry-polished with a 0.1 /im 
diamond micropolisher. The surface was black and shiny 
with an area of about 8 mm^. Because of ellipsomet- 
rically detected changes of the optical properties in air, 
probably due to surface oxidation, the sample was in- 
serted into the cryostat within minutes after polishing 
to make sure that the near-surface region corresponds to 
the bulk properties. Normal- incidence reflectivity mea- 
surements in polarized light did not reveal any optical 
anisotropy implying no preferential grain orientation. 

We combined ellipsometry and reflection spectroscopy 
to determine the optical constants for 15 K< T < 300 K. 
The dielectric function e = ei+ie2 was measured directly 
using a commercial optical ellipsometer with a UHV cryo- 
stat in the range 6000 - 37000 cm^^ (0.45 - 4.6 eV) at 
the angle of incidence 9 = 80° . The measurements at 6* = 
60° gave almost the same dielectric function, which is also 
an indication of isotropic properties of our sample. The 
normal incidence reflectivity spectra were measured from 
20 to 6000 cm-i (2.5 meV - 0.45 eV) utihzing Bruker 
113v FT-IR spectrometer, equipped with a home-design 
ultra-stable optical cryostat. Absolute reflectivities for 
each temperature were calibrated with the aid of a gold 
film, deposited in situ on the sample. In this range the 
complex reflectivity phase was obtained with Kramers- 
Kronig (KK) transformation, allowing us to calculate ei 
and 62. The phase spectrum was anchored using ellipso- 
metric data at higher frequencies (see Ref. p4[). 
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FIG. 2. Experimental spectra at 50 K and 300 K: normal 
reflectivity (left) and dielectric function (right). The left in- 
set shows the reflectivity in the FIR range for T =300, 250, 
200, 150, 100, 50 and 15 K (from bottom to top). The 
right inset shows the reflectivity in the full frequency range 
(at high frequencies, restored from ellipsometry data). The 
Drude-Lorentz fit at 300 K is shown in gray. 

An important issue is the effect of impurity phases on 
the optical spectra. The main impurity component is 
MgO, which has almost zero optical conductivity in the 
studied frequency range except for a few sharp phonon 
peaks in the FIR range. Using the Maxwell-Garnett the- 
ory of the two-component media p5| one can show that 
the role of impurity in this case is to reduce the ob- 
served conductivity uniformly by 7 - 10 % compared to 
the impurity-free sample. As impurity concentration is 
not precisely known, we do not correct the conductivity 
for this effect, but keep it in mind in the discussion. 
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FIG. 3. Optical conductivity at 50 K and 300 K determined 
from the experimental data (as described in the text). The 
left inset shows the T-dependence of the narrow Drude peak 
is shown. The right inset: the effective number of carriers per 
cell, N^g{u;) (solid curve: 50 K, dashed curve: 300 K). 



The experimental spectra are shown in Fig.y. At 
high energies the spectra demonstrate very little 
T— dependence. In the IR range the reflectivity steadily 
increases with cooling down. An additional small 
superconductivity-related increase of R is observed be- 
low Tc for w < 12 meV, in agreement with earlier data 
|16| ] . The main low- frequency features of the optical con- 



ductivity (Fig.H) are a narrow DP with a small spectral 
weight, and a broad "dome" structure extending up to 
1-1.5 eV. There is a minimum at « 0.15 eV and a 
broad maximum at « 0.3 eV. The structure itself looks 
like a strongly damped second DP. Another broad peak 
is centered at about 3.3 - 3.4 eV. One can note some 
phonon-like structures in the FIR reflectivity and, cor- 
respondingly, conductivity curves. Those structures are 
not reproducible from sample to sample, therefore we re- 
late most of them to impurity phases Ijl^l • 

The standard procedure of quantifying optical spectra 
is the Drude-Lorentz data fitting: e(w) = eoo — ^p/{oj{uj + 
ij)) + J2i ^pi/i^M ~ 1^^ ~ hi^), where flp and 7 are the 
plasma frequency and the scattering rate for the DP, and 
woi, ^pi and 7i are the frequency, the "plasma" frequency, 
and the scattering rate of the i-th Lorentz oscillator (LO), 
Coo is the high-frequency dielectric constant. We fitted 
simultaneously R{uj) at low frequencies and both ei(a;) 
and e2{^) at high frequencies (see Fig. ||) with one set of 
parameters listed in the Table for 300 K. 
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FIG. 4. Contributions to the optical conductivity in the 
Drude-Lorentz fitting at 300 K. The solid line is the experi- 
ment, the thick gray line is the fit; the dashed (dotted) lines 
show the narrow (broad) Drude peak and all the other curves 
show different LOs. The insets: temperature dependence of 
the Up and 7 for the narrow Drude peak. 

A good fit can be obtained by introducing two DPs 
with very different properties (Fig.H). The narrow 
DP has an almost temperature-independent plasma fre- 
quency « 1.4 eV; the temperature dependence of the 
scattering rate (see inset in Fig.^ scales well with the 
pDciT) curve (Fig|^(a)). The second DP is anomalously 
broad (7 w 1.1 - 1.2 eV). Although its plasma frequency 
(« 4.9 eV) is much larger than that of the narrow peak, 
the latter determines the major part of the static con- 
ductivity. The "dome" structure can be fitted with a 
combination of a broad DP and two LO at ~ 0.3 and ^ 
0.8 eV. Finally, there is another broad peak centered at 
^ 3.4 eV and a peak at ~ 5.1 eV. 

Let us now compare the data with the first-principles 
calculations |p|,p| . To account for the polycrystallinc sam- 



ple texture we averaged the in-plane and out-of-plane 
contributions to the optical conductivity using the effec- 
tive medium approximation (EMA) |l5[ assuming spher- 
ical randomly oriented grains. This approach is more 
accurate at low frequencies, where the wavelength A is 
larger than the typical grain size d (in our case of the 
order of 1-5 /im). For the limit A 3> <i the reflectiv- 
ity should be averaged directly rather than conductivity. 
For highly anisotropic compounds, e.g. cuprates, it was 
shown that some artificial structures may show up due to 
reflectivity-averaging effect 18 1. Although we take such 
a possibility very seriously, the fact that the complex di- 
electric functions obtained by ellipsometry measurements 
at incidence angles 60° and 80° match very well is an 
argument that the concept of effective medium conduc- 
tivity works sufficiently well even in the high-energy part 
of our measured spectral range. 

The interband optical conductivity (Figjs]) was calcu- 
lated by the WIEN97 LAPW code Q using the band 
structure of Ref. [^. The experimental peak at 0.3 eV 
clearly corresponds to the in-plane interband transition 
between the two cr-bands at 0.35 eV. The broad 3.4 eV 
peak finds no counterpart in the theory: the closest peak 
(2.4 eV, in-plane polarization) stems from the van Hove 
singularity at the M point. The reason for this discrep- 
ancy is not clear at the moment, although it might be 
related to the disorder effects discussed below. The in- 
crease of the conductivity above 4 eV, which is simulated 
by the LO at 5.1 eV, one can probably associate with the 
intensive out-of-plane maximum at 4.7 eV due to charge 
excitations between the B and Mg planes p9 . 

The intraband conductivity is particularly interesting. 
The calculated in-plane (out-of-plane) plasma frequencies 
of the 2D cr-band are 4.1 (0.7) eV, which in the EMA 
gives 3.0 eV for the effective plasma frequency uipa of the 
cr-band. Similarly, for the 7r-band, we obtain Wpjr = 6.2 
eV. The total plasma frequency tOp is thus expected to 
be about 6.9 eV. 

The observed flp of the narrow Drude peak alone (1.4 
eV), in agreement with Refs. |0,|[, is much smaller than 
the calculated total LOp. It can account for only a small 
part of the expected total spectral weight. However, 
about 80% of the missing spectral weight can be found 
in the features forming the broad "dome" structure. We 
suggest therefore that the narrow DP is due to the 2D cr- 
band carriers, while the broad "dome" is mainly formed 
by the intraband excitations within the 3D 7r-band with 
some contribution from the low-energy interband tran- 
sitions between cr-bands. The two orders of magnitude 
difference in the elastic scattering rates for the a- and 
TT-bands can be explained by the principally different na- 
ture of the two band systems. The cr-bands have hardly 
any weight on Mg, while the 7r-bands hybridize with Mg 
orbitals. On the other hand, the strong covalent bonding 
in the B planes makes any defect in this plane energeti- 
cally very unfavorable, while the loosely bound Mg planes 



should be prone to all kinds of defects. This implies that 
electrons in the cr-bands may experience much less defect 
scattering than those in the 7r-bands. 

Even this model cannot fully explain the small plasma 
frequency (1.4 eV) of the narrow DP, since the calculated 
cr-band effective plasma frequency is 3 eV. However, an 
important piece of physics is the mass renormalization 
due to strong EPI. This effect is neghgible if 7 » ujpk, 
which is usually the case in metal optics. In MgB2, how- 
ever, the EPI is believed to come mainly from the Fi2g 
phonon with ujph ^ 600 cm~^, larger than the a-band 
scattering rate. In this case, the mass renormalization 
still does not change the DC resistivity, but the peak 
gets narrower, and fitting by the standard Drude formula 
produces an cUp, which is, compared to the bare plasma 
frequency, reduced by a factor m* /m = 1 + A, where A is 
the EPI constant [gl|. The remaining spectral weight is 
transferred to higher energies, uj > cuph- The calculated 
couphng constants [|| for different bands are: Aqo- = 1.2 
and Aott = 0.45. This renormalizes the plasma frequency 
of the cr-band DP to « 2 cV. 

Given such a large scattering rate, is the Drude model 
applicable at all to the 7r-bands? The calculated Fermi 
velocity is about 4.9 x 10^ cm/s [|[ so that the 7 =1.1 
eV =2.7 X 10^^ s~^, results in a mean free path I = vp/'y 
~ 18 A, comparable with the lattice constant. This sig- 
nals approaching weak localization of charge carriers by 
defects, as opposed to Boltzmann transport. The grain 
boundaries might introduce additional scattering, how- 
ever, such a short mean free path is difficult to explain 
by the grain effects only. Localization effects should lead 
to deviations from the pure Drude formula: an additional 
suppression of conductivity at w — > and a characteris- 
tic maximum on cti (lo) curve at frequencies of the order 
of the inelastic scattering rate. This might explain the 
broad peak at 0.8 eV, which has almost the same 7 as 
the broad DP, and comparable spectral weight (see the 
Table) . The combination of these two terms might mimic 
the deformation of the conventional DP due to localiza- 
tion effects. 

Fig.|| shows the effective conductivity (including both 
interband and intraband components), obtained from 
first-principle calculations. We calculated the intraband 
absorption separately for the a— and tt— bands, and in- 
cluded the frequency-dependent mass renormalizations 
and scattering rates due to EPI, as described in Ref. |2^. 
The EPI Eliashberg function was taken in the Einstein 
form: a'^^F{uj) — 0.5Xoa-,TT^S{Lu — LOph)- The values for 
the T = scattering rates were taken from the experi- 
ment: 7ct= 15 meV, y^c = 1-1 eV. 

There is a qualitative agreement between the first- 
principles calculations and the experiment, despite a 50 
% difference in plasma frequencies. However, without 
the two-band scenario the discrepancy would be an or- 
der of magnitude bigger. The remaining error may be 
due to one or more of the following: above mentioned 7 



- 10 % scaling of the conductivity due to finite impurity 
concentration; surface quality issues; grain boundary and 
localization effects; possible underestimation of A in the 
calculations because of nonlinear EPI. 

Let us now make a link to superconductivity. Vari- 
ous experiments indicate that MgB2 cannot be described 
by the standard BCS theory |^. The popular two-gap 
theory p3,pl seems to be able to explain the majority of 
the experiments. However, it is known [|l3| that in this 
case the interband impurity scattering must rapidly sup- 
press Tc by isotropization of the order parameter. Such 
a suppression has not been observed [g|, which can be 
interpreted as a serious argument against the two-gap 
model. Our two-Drude picture of the low-energy optical 
conductivity helps to resolve this paradox. Indeed, the 
existence of two distinguishable Drude peaks would only 
be possible provided that the a — tt charge scattering is 
small. In this case impurity scattering, albeit reducing 
the DC conductivity, does not isotropize the gap, there- 
fore the Anderson theorem holds and no pairbreaking 
occurs Q. 

It is worth noting that one cannot, based just on the fit 
quality, choose between our "two-DP" model and a DL 
fit with the second DP replaced by a broad low-energy 
(■^0 ^ 0.25 eV) LO. Although the latter interpretation is 
consistent with the data, we favor the two DP model for 
the following reasons: (i) the plasma frequency, unlike 
some other electronic properties, is known to be repro- 
ducible by LDA calculations within errors hardly ever 
larger than 100 %, as opposed to a factor of five discrep- 
ancy we would have for MgB2 . (ii) With one DP the total 
u)p is less than 2 eV; in this case the application of the 
standard high-T limit of the Bloch-Griineisen formula, 
P'Dc{T) = Sn'^XtrT/i^p would give the EPI transport con- 
stant Xtr < 0.15 g, an unphysically small value for an 
covalent-ionic metal with a sizeable density of states. Fi- 
nally, (iii) the two-DP model resolves naturally the main 
problem of the two-gap superconductivity model, sup- 
ported by a variety of experiments. 
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FIG. 5. Comparison of the experimental and calculated op- 
tical conductivity, including the interband part shown in the 
inset, at 300 K (see text). 



In summary, the optical conductivity of MgB2 in a (1995). 

wide frequency range was obtained using a combination 

of ellipsomctry and normal incidence reflectomctry meth- 

ods. The drastic suppression of the low-energy spec- ^ABLE I. Parameters (in eV) of Drude (1, 2) and Lorentz 

, 1 • 1 , c J i 1 1 i- ■ J.^ ij.-i_ J (3-6) peaks obtained by spectra fitting at 300 K (eoo = 2.9). 

tral weight fands a natural explanation m the multiband y ' t- j t- & y oo y 

model: a narrow Drude peak corresponds to the 2D a- 
bands (additionally renormalized by the EPI) , while the 
intraband absorption of the 3D 7r-bands is spread over a 
large frequency range up to 1 - 1.5 cV. Optical measure- 
ments on high-quality single crystals of MgB2 are neces- 
sary to get more reliable data on the optical conductivity 
and its anisotropy. 
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